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Many topological phenomena first proposed and observed in the context of electrons in solids 
have recently found counterparts in photonic and acoustic systems. In this work, we demonstrate 
that non-Abelian Berry phases can arise when coherent states of light are injected into “topological 
guided modes” in specially-fabricated photonic waveguide arrays. These modes are photonic ana¬ 
logues of topological zero modes in electronic systems. Light traveling inside spatially well-separated 
topological guided modes can be braided, leading to the accumulation of non-Abelian phases, which 
depend on the order in which the guided beams are wound around one another. Notably, these 
effects survive the limit of large photon occupation, and can thus also be understood as wave phe¬ 
nomena arising directly from Maxwell’s equations, without resorting to the quantization of light. 
We propose an optical interference experiment as a direct probe of this non-Abelian braiding of 
light. 


Manifestations of topology in physical systems, specif¬ 
ically in the form of so-called geometric phases, mm 
have risen to prominence over the last three decades. 
Geometric phases were shown by M.V. Berry to arise 
in quantum systems under cyclic adiabatic variation of 
parameters, [2] but it was later understood [3] that this 
phase had been discovered thirty years earlier in the 
context of classical optics by S. Pancharatnam. [1] The 
close analogy between quantum mechanics and classi¬ 
cal optics has remained over the years, and many of 
the striking topological states of matter associated with 
electrons in solids, such as topological insulators and 
semimetals, [4] have recently found counterparts in pho¬ 
tonic [SHIS] and acoustic inHis] systems. The present 
work aims to further highlight and deepen the connec¬ 
tion between topological phenomena in solids and in 
classical wave mechanics by demonstrating a new facet 
to this correspondence. We demonstrate the existence of 
a classical analogue of the non-Abelian Berry phase [20] 
that arises from “braiding” topological defects in solids. 

One route to non-Abelian Berry phases in electronic 
systems lies in the remarkable physics of zero modes. 
As was pointed out in pioneering work by R. Jackiw 
with C. Rebbi [21] and P. Rossi, [22] localized zero- 
energy fermionic states can bind to topological defects 
in an order parameter, such as kinks in one dimen¬ 
sion and vortices in two dimensions. In systems where 
electric charge is a good quantum number, these zero 
modes carry charges that are fractions of the “funda¬ 
mental” electron charge. Eu usi m In chiral super¬ 
conductors, where charge conservation is broken, these 
localized modes are Majorana bound states with non- 
Abelian statistics. [22] [251127] At the mean-field level, 
where interactions between electrons are neglected, this 
non-Abelian statistcs can be understood in terms of 
the accumulation of non-Abelian Berry phases as de¬ 
fects are adiabatically exchanged with one another. [28] 
The non-Abelian nature of this process manifests itself 
in dependence of these phases on the order in which 


the defects are interchanged, in contrast to the usual 
“Abelian” Berry phase. 

In this Letter, we propose a novel means of realizing 
topological zero modes in photonic, rather than elec¬ 
tronic, systems, and demonstrating their non-Abelian 
braiding directly. The proposed realization consists of 
non-interacting photons propagating in the non-trivial 
background of a photonic lattice with topological de¬ 
fects whose positions are controllable. Light channeled 
into “topological guided modes” localized at these de¬ 
fects can be braided, leading to the accumulation of non- 
Abelian phases that depend on the order in which the 
braids occur. We demonstrate that this effect manifests 
itself at both the quantum and classical levels, owing to 
the linearity of the equations of motion for noninteract¬ 
ing light. 

The zero modes we propose to realize are photonic 
analogues of Kekule zero modes in graphene, [24] which 
are bound to vortices in the complex order parameter 
A(r) describing a dimerization pattern in the hexagonal 
lattice (Fig. I). The translation between electrons and 
photons is achieved by replacing the sites of the lattice 
with waveguides embedded in a bulk optical medium 
{e.g. fused silica), [29] which are extended in the 2 : di¬ 
rection and whose x-y positions mimic the 2D posi¬ 
tions of the carbon atoms in the distorted graphene lat¬ 
tice. The wave equation for the paraxial propagation 
of light in such a waveguide array maps directly onto 
the time-dependent Schrodinger equation (SE), where 
the time coordinate t in the SE is replaced by the coor¬ 
dinate ^ along the direction of light propagation. This 
wave equation can be further mapped, using coupled 
mode theory, m to a linear differential equation that 
is in one-to-one correspondence with the noninteract¬ 
ing tight-binding model of the electronic system. The 
waveguides themselves can therefore be thought of as 
the world lines of the carbon atoms, with straight waveg¬ 
uides corresponding to a lattice that is static in time. 

Waveguide arrays of this type have been realized. 
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FIG. 1: a) Uniform Kekule distortion in a hexagonal waveg¬ 
uide lattice at a slice of constant z. The faint circles repre¬ 
sent the original hexagonal lattice, while darker circles repre¬ 
sent the distorted lattice. Sublattice A is colored red, while 
sublattice B is colored blue. The overlaid vector field rep¬ 
resents the magnitude and direction of the order parameter 
A(r). The background is an image of a low-lying eigenmode 
of Eq. Q . The intensity pattern represents the amplitude of 
the electric field in each waveguide, b) Waveguide lattice in 
the presence of a vortex in the Kekule pattern, with an over¬ 
laid order-parameter vector field showing circulation around 
the vortex core. The background is an intensity plot of the 
localized zero-mode wavefunction. We have chosen a vortex 
profile in which only sublattice A is displaced. 


written into bulk materials with exquisite precision us¬ 
ing femtosecond lasers, m (see Ref. |29| for a review). 
They have been used to mimic graphene and other elec¬ 
tronic systems with both static lattices O I32H34] and 
ones that change as a function of time. [TOj [11] The 
experimental protocol suggested here hinges on this ca¬ 
pability to execute a braiding procedure in which three 
vortices are wound around one another. As we show 
below, a protocol can be chosen that reveals the non- 
Abelian nature of the braiding directly via interference 
between different braiding patterns. 

Our starting point is the coupled-mode equation for 
paraxial light propagation through the waveguide array, 
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i dz ~ ^ ^ (f) 

i=i 

Here, the vector r is defined on a hexagonal lattice di¬ 
vided into two triangular sublattices that we call A and 
B. The vectors -\-Sj (j=l,2,3) connect the point r in 
sublattice A to its three nearest neighbors in sublattice 
H, located a distance a away, and —Sj connect points in 
B to their neighbors in A. For simplicity, we assume the 
light in the array to be monochromatic, so that Eq. 0 is 
frequency-independent. Eq. 0 is formally identical to 
a noninteracting tight-binding model in which the time 
coordinate t has been replaced by the longitudinal co¬ 
ordinate 2 ;. Its solutions can be decomposed into eigen- 
modes ^ where = 0,..., 2N — 1, 

with N the number of unit cells in the waveguide array, 
labels the “energy eigenvalues” The “wavefunction” 


ipr{^) ill Eq. 0 is related to the amplitude of the local 
electric field Er{z^ t) on the “site” r. It is useful to work 
with the quantized electric field operator [35] 

Er.{z, t) = Y^ b^n + H.c. (2) 

V 

where n is a unit vector in the x-y plane describing the 
polarization of the field and satisfies the dispersion 
relation of a single waveguide. The ladder operators 
bjy satisfy the bosonic commutation relations = 

[See Supplemental Material (SM) for details on 
the quantization procedure.] 

The “Hamiltonian” in Eq. 0 depends exponen¬ 
tially on distances between nearest-neighbor waveg¬ 
uides (we neglect longer-range couplings for the mo¬ 
ment). We take Hr^r+sj = —t — Strj, where 6trj = 
A(r) -|- c.c. (we denote by K± the lo¬ 

cations of the two inequivalent Dirac points, at oppo¬ 
site corners of the hexagonal Brillouin zone). Here, the 
parameter t describes the evanescent couplings of the 
waveguides, and the position-dependent function Strj 
describes modulations of these couplings due to dis¬ 
placements of the waveguides from their original x-y 
positions. The complex-valued Kekule order parame¬ 
ter A(r) controls the distortion of the lattice. [24] A 
vortex in A(r) is a defect in this distortion pattern, 
but not the lattice itself. The order parameter in the 
presence of a single vortex centered at the origin reads 
A(r) = Ao{r) in polar coordinates r = {r,0). 

Here, Ao(r) = A tanh(r/^o) describes the spatial pro¬ 
file of the vortex, which has a core radius and a is 
the phase of the order parameter. 

The zero mode in the presence of this vortex profile 
can be found by setting the left-hand side of Eq. 0 
to zero and solving for 'ijjr. The zero-mode solution is 
tightly localized near the core of the vortex, with a size 
of order 1/A, and has support on sublattice A only [see 
Fig. [^B)]. [If we send 0 ^ —0 in A(r), the zero mode 
has support on sublattice B instead. [24]] This means 
that light propagating in the zero mode travels as if con¬ 
fined to an “optical fiber” located at the vortex core, al¬ 
beit with evanescent decay into neighboring waveguides 
in the same sublattice. However, the zero mode differs 
crucially from a mode in an optical fiber, both because 
it takes the distortion of an entire waveguide lattice to 
create, and because it depends on the topological nature 
of this distortion. These “topological guided modes” are 
responsible for the non-Abelian effects described below. 

The above discussion neglects the presence of various 
lattice effects that can shift the zero mode eigenvalue to 
a finite momentum. In principle, this makes possible the 
accumulation of path-dependent {i.e. non-topological) 
dynamical phases during braiding. However, as we show 
in detail in the SM, these effects can be neglected to 
resonable precision without fine-tuning. There, we es¬ 
timate (conservatively) that we can perform thousands 
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of braids before the effects of dynaniical phases begin to 
manifest themselves. 

To facilitate our discussion of the non-Abelian braid¬ 
ing of vortices like the one above, we work in the con¬ 
tinuum limit, where the analytical form of the vortex 
wavefunction is known. (Exact diagonalization calcu¬ 
lations corroborating these results are reviewed in the 
SM.) On length scales much longer than a, the operator 
that annihilates a photon in the zero mode is 


^0 = y dr 


u{r) 6 , (r) -h b_{r) 


(3) 


where b±{r) annihilate a photon at the Dirac points 
K±. (These operators are assumed to be normalized 
such that the canonical commutation relation = 

1 holds.) The function u{r) = u{r) is given, up to nor¬ 
malization, by m 


and the zero-mode Hilbert space is spanned by the op¬ 
erators ^ A clockwise adiabatic exchange of vortices 
i and i + 1 is implemented by winding the vortex-core 
coordinates Ri{z) and Ri-^i{z) around one another as 
functions of increasing 2 ; (see, e.g., the braids in Fig. |^. 
In the SM, we show that the nontrivial nature of this 
winding process stems from the double-valuedness of 
u{r) under a ^ a — 27r. (We verified this statement 
for the zero-mode wavefunction on the lattice via a nu¬ 
merical tight-binding calculation, as we describe in the 
SM.) The effect of this exchange, up to a gauge choice, 

is to map b^- ^0 i+i ^0 i+i ^ ~^o i while leav¬ 

ing all other vortex operators unchanged, similarly to 
the case of Majorana [26] and Dirac zero modes [36] in 
electronic systems. The operator that implements this 
exchange is 






( 8 ) 


^ gi(a/2+7r/4) dr' Ao(r') _ (^4^ 

Note that the zero mode amplitude u{r) is double¬ 
valued: when (a, the phase of the order parameter A(r), 
changes by 27r, u{r) acquires a minus sign. The origin 
of this double-valuedness lies in the fact that u{r) is 
a solution to the Dirac equation, and therefore must 
transform as a spinor under changes in a. 

The quantum states created by 6 q can be connected 
to the macroscopic state of light in the waveguide array 
by defining the coherent state 

|^) = g-|A|V2gA6t|o) 

with mean photon number (SoSg) = |Ap. These states 
provide a faithful description of the electric field in the 
waveguide array in the large-|A| limit when the input 
light is a coherent, monochromatic laser spot centered 
on the vortex core. In this case, the contribution of the 
photonic zero mode to the classical electric field is given 
by [c.f. Eq. ([ 2 ])] 

EoAz,t) = AM(r)e*(*^+-’’+'=-^-“‘)n + c.c., (6) 

with u{r) given in Eq. Q and A defined by Eq. ^ . The 
ease with which one can translate between the quantum 
and classical descriptions of the waveguide array owes 
to the fact that photons in the array are noninteracting. 
Indeed, an alternative way of deriving Eq. (§ is to solve 
Maxwell’s equations directly in the presence of a vortex. 

Let us now study the braiding of zero modes in an 
infinite system with v vortices at z- dependent positions 
Ri{z) (i = 1,..., 1 ;). The order parameter in the pres¬ 
ence of this vortex configuration is given by 

n 

A(r) = A tanh(|r - Rj\/£o) , (7) 

i=i 


One verifies by direct calculation that these operators 
satisfy = 0 for \i — j\ > 1 and TiTjTi = TjTiTj for 

— jI = 1, and therefore form a unitary representation 
of the braid group on v strands. m The action of 
these generators on a state 


|ni,.. 


)=n 


(biiT 




| 0 ), 


(9) 


representing a system with v vortices, each with a fixed 
number of photons, is given by 


fi\.. ... ) = (-l)”’+i|... ,ni+i,ni,...). (10) 


We now explain the consequences of the above phase 
factor, which arises purely from braiding the vortices. 
Consider the case of two vortices, and the operation of 
winding one around the other, which restores them to 
their initial locations. We begin by considering i) an 
example where this phase factor does not lead to uni¬ 
tary operations on a multi-dimensional space of degen¬ 
erate states, and then turn to ii) an example where it 
does. For case i), consider any eigenstate of occupation 
number, |ni,n 2 ), which upon winding goes to the state 
(^_l^ni-\-n 2 \ni^n 2 ); clearly, the initial and final states are 
equal up to a phase, and the braiding operation acts on 
a one-dimensional space only. In contrast, consider case 
ii), where we take an initial state that is not an eigen¬ 
state of occupation number, but rather a superposition 
of states with occupation 0 and 1, for example. Then 
winding takes 


|0) + |1) |0) + |1) , |0)-|1) |0)-|1) 
^/2 ® V 2 ^/2 ® ,/2 


( 11 ) 


Clearly, the initial and final states are different (here, 
they are orthogonal). This is a crucial hallmark of the 
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non-Abelian nature of the action of the braiding gen¬ 
erators on these states. Thus, the essential ingredi¬ 
ent necessary for braiding to connect states in a multi¬ 
dimensional Hilbert space is that the initial state of the 
zero mode system consists of superpositions of states 
with even and odd numbers of photons. 

Such superpositions do not only exist within the do¬ 
main of quantum optics. The coherent states defined in 
Eq. (§ consist of a superposition of photon states with 
all occupation numbers, and can be created by shining 
a coherent laser beam centered on a single vortex core. 
Let us now consider a system of v vortices, into each of 
which is loaded a coherent state of photons: 

= ( 12 ) 


The action of the braiding generators on these states is 
found from Eq. to be 


• • • : ^25 • • • ) — I • • • 5 ^25 • • • )• (13) 


To see that this braiding operation connects quantum 
states in a multi-dimensional space of degenerate states, 
consider again the winding of two vortices, which re¬ 
stores their initial locations. Eor a coherent state, this 
operation takes 


1^1, A 2 ) —^ I — Ai, —A 2 ) . (14) 

The overlap of these two states is (Ai, A 2 I — Ai, — A 2 ) = 
g- 2 (|Ai| +IA 2 I ) ^ whose magnitude is smaller than one, in¬ 
dicating that the states are linearly independent. This 
demonstrates that the Hilbert space spanned by the 
braiding operations is multi-dimensional. Moreover, in 
the limit where a large number of photons are loaded 
into the zero modes (large |Ai^ 2 |: which is to be expected 
from a laser), the overlap is exponentially small, and the 
initial and final states become orthogonal. 

At this point, it is worthwhile to reflect on the signif¬ 
icance of the fact that the effect of braiding the vortices 
manifests itself at the level of coherent states, which 
are essentially classical. While the braiding of these 
vortices can be understood at the quantum {i.e. few- 
photon) level, its effects permeate the entire zero-mode 
Hilbert space for arbitrary occupation numbers. Conse¬ 
quently, the quantum action of the braiding generators 
survives the limit of large occupations, so that macro¬ 
scopic effects of this braiding can be seen. In particular, 
observe that Xi —Xi under a braid corresponds, in the 
limit of large | A^ |, to a change in the sign of the electric 
field Ei near the core of the ith vortex [c.f. Eq. §]• This 
observation forms the basis of our discussion below. 

We now propose an experiment that could provide di¬ 
rect and unambiguous evidence for non-Abelian braid¬ 
ing at the level of coherent states. The photonic non- 
Abelian interferometer (PNAI) that we propose, shown 



FIG. 2: Schematic of the proposed photonic non-Abelian 
interferometer. A coherent laser beam (one of three, one 
for each vortex in a lattice) passes through a 50/50 beam 
splitter (blue) and simultaneously enters two separate pho¬ 
tonic lattices fabricated with the two braids to be compared. 
The two output beams are reflected by mirrors (grey) into 
another beam splitter that interferes the two signals and out¬ 
puts the sum and difference to separate screens (black). For 
the choice of braids used here, one screen shows a bright 
and two dark spots, while the other screen shows the “logi¬ 
cal complement” of the first, with one dark and two bright 
spots. 


in Fig. consists of two separate photonic lattices, 
each with three vortices, with waveguides written into 
the host medium in such a way that the vortex cores 
wind adiabatically around one another according to a 
specific pattern. (The adiabatic condition here corre¬ 
sponds to demanding that any change in the position of 
each waveguide as a function of 2 : occur on length scales 
much larger than the inverse photonic bandgap.) In one 
lattice, the waveguide pattern executes a braid H 1 H 2 , 
while in the other, the waveguide pattern implements 
the braid B 2 B 1 . Interfering the light output from the 
two lattices reveals the defining feature of non-Abelian 
braiding, namely that performing the same braids in 
different orders yields different results. We now proceed 
through each stage of the interferometer setup. 

We start with the input stage of the PNAI. We assume 
that the input light comes from three monochromatic, 
coherent laser sources, each focused on the core of a 
single vortex so that there is large overlap with the zero 
modes. Each of the three input beams is split by a 
50/50 dielectric beam splitter, so that the light entering 
each vortex core comes from the same source beam. We 
will denote the light entering the upper branch of the 
interferometer in Eig. [^as |Ai,A 2 ,A 3 ), and |Ai,A 2 ,A 3 ) 
for the lower branch, ^he first beam splitter enforces 
the phase relation Xj = —iXj.) 

Eor the braiding stage of the interferometer, we 
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choose, for example, Bi = T 2 T 1 and B 2 = T 2 T^^T 2 T^^ 
(see Fig. [^, which ensures that the vortices on the out¬ 
put facets of both lattices are in the same order. [38] 
In the SM, we provide more information about how to 
write these braiding patterns into the waveguide array. 
The output states from the two braids are 

B1B2 |Ai, A2, A3) = I — Ai, A2, —A3) (15) 

and 

B2B1 |Ai, A2, A3) = I — Ai, —A2, A3). (16) 

In the final stage of the interferometer, the output 
beams from the braiding stage are combined at another 
beam splitter. The sign differences between the coher¬ 
ent states exiting the two branches of the interferometer 
cause the light to interfere constructively at one detec¬ 
tor and destructively at the other. Which detector this 
is for each vortex depends on the relative signs of ^ 
and Ai, which in turn depend on the braids (see Fig. ^ 
and SM for more details). Since the effects of dynamical 
phases can be heavily suppressed in a controlled manner, 
as discussed earlier, the only source of this interference 
is the noncommutativity of braiding the vortices. 


In summary, we have demonstrated in this Letter a 
means to realize photonic analogues of topological zero 
modes in photonic lattices. We demonstrated that these 
“topological guided modes” can be understood at both 
the quantum and classical levels when the photons in 
the waveguide array are weakly interacting. We fur¬ 
ther proposed a photonic non-Abelian interferometer, 
feasible with current technology, to detect unambiguous 
signatures of the non-Abelian Berry phases that result 
from braiding these topological guided modes. 
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Supplemental Material: Non-Abelian braiding of light 


Appendix A: Quantized electromagnetism in a waveguide lattice 


We seek to quantize the electromagnetic field in the presence of a spatially varying relative permittivity e(r), 
and follow the notation of Quantum Optics^ by Vogel and Welsch. [35] (We adopt the convention here that r is 
a three-dimensional vector when it enters the argument of a function, and a two-dimensional vector in the x-y 
plane when it appears as a subscript.) We denote the vector potential by A{r,t), and work in the Coulomb gauge 
V • A = 0. In the absence of free charges we also have Aq = 0. We work under the assumption that Ve(r) • A is 
small. With this the classical Hamiltonian is 

H = [d^r(+—[V X A]A (A1) 

J \e{r)eo jio J 

with n(r,t) the momentum conjugate to A{r^t): 

n(r, t) = e(r)eo A(r, t). (A2) 

We then have the canonical Poisson brackets: 

{A{r,t),U{r',t)} = S{r-r'). (A3) 

We quantize Eq. ( |A1| ) by promoting A and 11 to operators obeying equal-time commutation relations 

[A(r, t), n(r', t)] = ihS{r — r') (A4) 

with all other commutators vanishing. Our operator equations of motion, 

A(r,t) = ^[A{r,t),H(t)\ = ( (A5) 

in e(r)eo 


in jiiQ 


lead to the second-order equation of motion for A(r^t) : 


V^A(r,t) 


e(r) d‘^A{r, t) 


(A6) 


(A7) 


For a quasi-two-dimensional array of waveguides, our relative permittivity takes the form e(x, z) = e^[l+Ae(x, y)]. 
We can solve this using separation of variables with A{r^t) = ^ Ajy^^{r)Tjy^^{t)bjy^^ -j- H.c., where we are 

anticipating two indices u and uj which will label solutions of the resulting equations for A(r) and T{t): 


Tu,uj{t) + oj'^ = 0 


(A8) 
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and 


^9 ^ / N oo‘^er ^ , X uj‘^erAe(x,y) ^ . 

\/^Ai,^u{x,y,z) H-H-^- —A^^^{r) = 0. 


(A9) 


The first equation is easily solved with = T(0)e“*^^, and restricting cj > 0. The second equation is the 

familiar Helmholtz equation, which the classical vector potential would also obey. We are interested in waves 
propagating primarily in the positive ^-direction, and thus consider solutions of the form 

y, z) = -iuj~^ y, z) n. (AlO) 

Here v = cj and the factor has been chosen so that here corresponds to in equation (2) in the 

main text. This leads to an equation of motion which is first-order in 2 ; (assuming the envelope varies slowly), 


idz^u,uj{x,y,z) 


1 

2 



-1 


1 

2 



Ae(x,y) ^^^^(x,y,z), 


(All) 


where is the Laplacian operator acting in the x-y plane. This resembles a two-dimensional Schrodinger equation 
with the ^-direction taking the place of time, and with a potential Ae(x, y). Motivated by this analogy, we translate 
the Schrodinger-like equation into a lattice model by writing ^{x^y^z) (we will suppress the u and uj indices for 
now) as a sum over “Wannier states” of the lattice: ^{x^y^z) = '^^'ipr{z)wr{x,y). Using the orthonormality of 
Wannier states we find 


idz'tpriz) = 'Ipriz) {Wr\H\Wr) + '^r'(z) {Wr\H\Wr'), 

r'^r 


(A12) 


where we have switched to bra-ket notation for the Wannier states, and H = —^(^)~^V^ — |(^)Ae(x,^). Matrix 
elements between different Wannier states decay exponentially with the spacing between lattice sites, and so we 
neglect all off-diagonal elements except Hr^r+sj = {xOr\H\wr+sj)• We remove the diagonal elements, which adjust 
the wavevector to satisfy the massive dispersion relatio n of t he waveguide at the “site” r, by defining = 
uj/v + {wr\H\wr) aud replacing uj/v in equation ( |A1Q[ ). With these adjustments, we arrive at the tight- 

binding equation for 'ipriz), 


idz'lpri^Z^ — ^ ^ (Al3) 

j = l,2,3 

Here the sign choice is positive if r is on sublattice A and negative if on sublattice B. Restoring indices, and now 
using u to index eigenvectors of this Hamiltonian, we have 

y,z) = '^ Wr{x, y). (A14) 

r 

where ^u,uj,r is an eigenstate of the Hamiltonian H^ r,r±sj with eigenvalue Our full expressions for the vector 
potential and electric field are 

A(r, t) = ^(-icc;“^) ^ Wr{x, y) n + H.c., (A15) 

u,uj r 


E{r, t) = = y y Wr{x, y) n + H.c. (A16) 

With these solutions in hand, we can derive commutation relations of the and bl ^ by enforcing our commutation 
relations for A and 11. As expected, we find that: 


[^iy,uj ,uj'] — 0 


(A17) 




(A18) 
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(A19) 

Lastly, it is convenient to decompose E{rA) in the Wannier basis, and define an 
lattice site using: 

electric field operator at each 

^(^: t) Wr{x, y), 

r 

(A 20 ) 

Er{z, t) = Y^ „ + Yi.C. 

(A 21 ) 


iy,uj 


In the main text, we have specialized to the case of monochromatic input light, which selects a single frequency uj 
in the sum above. 

For the sake of simplicity, we have neglected in the main text the r-dependence of the waveguide dispersion 
This is rigorously justifiable when the variation in the relative permittivity Ae(x, y) has the discrete translational 
symmetry of the undistorted hexagonal lattice, or that of the translationally-invariant Kekule pattern shown in 
Fig. 1 (A). However, when translational symmetry is broken by the presence of a vortex, is not independent of 
the waveguide position r. This is problematic for systems with multiple vortices, as it can lead to the accumulation 
of dynamical phases as the vortices are braided (see Appendix E). We argue in Appendix E that the effects of these 
dynamical phases on the zero modes can be eliminated for an appropriate class of distortions of the waveguide 
lattice, presented in Appendix C. 


Appendix B: Zero-mode statistics and exact diagonalization results 

1. Hamiltonian 


We begin by reviewing the geometry of the hexagonal waveguide lattice. The vectors connecting a site on 
sublattice A with its three nearest-neighbors are 


Si = (0, -a), S2 = ^a), and S 3 = a, ^ a) 


(Bl) 


A unit cell consists of one site on sublattice A and a single nearest-neighbor on sublattice B. Our lattice vectors 
are then 


\/3 3 

ai = S 2 - Si = (— a, - a), 02 = S 2 - S 3 = (\/3a, O). 


(B2) 


We have chosen a to be the distance between nearest-neighbors, whereupon the distance between unit cells is '/S a. 
In reciprocal space, the Dirac points lie at opposite corners of the hexagonal Brillouin zone, K± = ± (, O). 
The Hamiltonian for the class of waveguide lattices that we study is given by 


H, 




— t Stj-' 


(B3) 


Hr-\-Sj ,r - Hrp (B4) 

for r in sublattice A and all other matrix elements zero. In the presence of a Kekule texture, 6tr,i is given in terms 
of the complex-valued order parameter A(r) by: 

5tr.,i = ^ A(r) e^(K+-K_)-r 

A vortex centered at position R is described by A{r;R) = Ao(|r — where Ao(|r — i^|) is a 

real-valued function which vanishes at \r — R\ =0 and approaches a constant value when \r — R\ ^ io, for some 
length scale Iq. Eor simplicity, we take Ao(|r — i^|) = A tanh(|r — i^l/^o), as in the main text. Eor multiple vortices, 
the order parameter is proportional to the product of order parameters for each individual vortex; for v identical 
vortices at positions Ri ,..., Ry , we have 


A(r) = A 


PJtanhdr - Rj\/io) 


i=i 


g* EJ=i -arg(r-Rj )] 


(B 6 ) 
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2. Vortex statistics 


The Kekule texture with v vortices admits v zero modes, one localized at each vortex. We assume vortices are 
kept far enough apart to prevent mixing betw een zero modes (see Appendix E). Near a single vortex, say near i^i, 
we can approximate r ^ Ri in equation (B6) for all contributions from vortices Rj^j > 1. So the zero mode at the 
first vortex sees an order parameter 


A(r) « Atanh(|r — Ri\/£o) 


]^tanh(|i^i — Rj\/io) 

j=2 


g-i aTg{r-Ri ) -arg(Ri-Rj)] 


Atanh(|r - R^\/£Q)e-^^'^s{r-R^)^^a^+iT.]=2l<^i-^rgiR^-R^)]^ 


(B7) 


where we have used that, for Iq much smaller than the distance between vortices, we have tanh(|i^i — Rj\/io) ^ 1 
for all i,jf. Thus, the only effect of the n — 1 other vortices on the first zero mode is to shift the phase of its order 
parameter by a constant, i.e. taking ^ o^i + ~ 

Now consider braiding two vortices at positions Ri and R2 in a counterclockwise fashion. The braiding operation 
takes arg(i^i — i^ 2 ) ^ arg(i^i — i^ 2 ) Ttt and arg(i ^2 — ^ 1 ) —> arg(i ^2 — ^ 1 ) Ttt, while leaving all other arg(i^i — Rj) 
unchanged. Using the local approximation of the order parameter near Ri in equation (B7), we see that braiding 
vortices Ri and R 2 is equivalent to taking ai ai — ir in the single vortex wavefunction. Using a similar 
approximation near the vortex at R 2 , we also have 0^2 ^ 0^2 — tt. Relabeling the exchanged vortices 1 2 and 

using arg(i^i — R 2 ) = arg(i ^2 — ^ 1 ) + tt, we finally have ai 0^2 and 0^2 ^ — 27r. Because the zero mode 

wavefunction (4) is double-valued with respect to a, we have u{r; a — 27r) = —u{r; a), and we arrive at the statistics 
in Eq. (8). 


3. Tight-binding braiding simulation 


We here review exact diagonalization calculations of the tight-binding Hamiltonian (B3)-(B6), to verify the 


existence of the zero modes and the action of the braiding generators in Eq. (8). Upon introduction of a constant 
A(r) = A, a band gap opens among the bulk states, although a small number of edge states with zero and near-zero 
energy remain. The introduction of a vortex in the order parameter leads to a single zero mode localized at the 
vortex core, and an additional edge mode, which are decoupled from one another and at zero energy for sufficiently 
large lattices. The zero mode has support only on sublattice A, and decays exponentially with distance from the 
vortex core with a characteristic length inversely related to |A|. 


Again exploiting the mapping between the Helmholtz equation (A9) and the tight-binding equation ( A13 ), we 
simulate the effect of braiding in our waveguide array by promoting our tight-binding Hamiltonian (B3) to a 
function of time, H{t) (recall that “time” in the tight-binding equation is actually the spatial 2 ;-direction in the 
waveguide array — we refer to it as t in this section to make it clear that our simulations were performed in a two- 
dimensional system), implicitly through its dependence on A{R) and thus the vortex centers Ri{t). We simulate 
a system where A(r) hosts two vortices, at positions Ri and R 2 . To perform a single counterclockwise (CCW) or 
clockwise (CW) braid of the vortices, we describe the vortex positions as a function of time by 


Ri{t) = Rb (cos(7rt/T), ±sin(7rt/T)), 


(B8) 


R 2 {t) = —Rb (cos(7rt/T), ±sin(7rt/T)), (B9) 

where Rb i^ the radius of the semicircle traced out by each vortex during the braid (or equivalently, half the distance 
between the vortices), T is the total time of the process, and the sign of the ^-component is determined by the 
handedness of the braid. The initial state of the system was taken to be an equal-amplitude superposition of the 
zero mode states centered at each of the vortices. We then discretize the time interval into N steps of duration 
dt = T/A, and act with the matrix on the state at each time tj = j 6t^ j = 1,..., A, achieving the final 

state after tj = N 6t = T. 

Results for a simulation on a 44 x 45-unit-cell lattice with Rb = 12 a, T = 120/t, A = 100, and A = .71 are 
shown in Eig. (§. The braiding generators (8) predict that under a CW braid, the vortex on the right undergoes 
a sign change between the initial and final states, while the vortex on the left does not. Eor a CCW braid, on 
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\i^i,R) \'(PcW,l) 


\'4^cw,r) 




b) 


0.2 


0.1 


0.0 


- 0.1 


- 0.2 

FIG. 3: (a) The superposition of vortex zero modes used as the initial state for both braids. Black dots are placed at the 
vortex centers, and a dashed line is drawn at the boundary used to decompose the initial and final wavefunctions into left and 
right pieces, as discussed in the text. To display the results in a more visually-friendly manner, the tight-binding wavefunction 
amplitude was multiplied by an arbitrarily chosen Gaussian “Wannier function” at each lattice site, (b) Results of the tight- 
binding simulation of a single GW braid. Arrows indicate the handedness of the braid performed. A visual comparison to 
(a) shows that the braid has resulted in a left vortex of the same sign as the initial state, and a right vortex of opposite sign, 
confirming the action of the braiding operators (8). (c) Results of the tight-binding simulation of a single GGW braid. The 
left vortex is now of opposite sign as the initial state, and the right vortex is of the same sign. 



the other hand, the vortex on the left undergoes a sign change, while the vortex on the right does not. A visual 
comparison of the final wavefunction after each braid to the initial wavefunction confirms this action (see Fig. |^. 
To quantify this, it is helpful to decompose each wavefunction as \ 2 p) = I^Pl) + \iPr), where has support 
only to the left of a vertical line equidistant from both vortices, and I'lpn) has support only to the right. We 
write the initial state of both braids as l^ipi) and the final state after CW/CCW braids as \'ipcw) and I'lpccw)^ 
respectively. In the ideal limit of well-separated vortices and adiabatic braiding, the braiding generators predict 
that {'ipi,L\'^cw,L) = {'^i,R\'^ccw,R) = 0.5, and ('0i,i?|'0cvp,i?) = {^i,L\^ccw,L) = -0.5. Note that the overlaps 
have magnitude less than unity because the I^jl) and I'iPr) are not normalized states. Computing the overlaps 
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numerically, we find 


{i^i,L\i^CW,L) = 0.47, 

(BlOa) 

= -0.39, 

(BlOb) 

('0i,L '0 ccw,l) = -0.38, 

(BlOc) 

{i^i,R\i^ccw,R) = 0.46, 

(BlOd) 

^^-10“^^) for all overlaps. Note that. 

as predicted, the vortex on the 


with zero imaginary part (to one part in 10 
left does not undergo a sign change under a CW exchange, while the vortex on the right does [c.f. Eqs. (BlOa) and 


(BlOb)]. Similarly, as predicted, the vo rtex o n the left ac quires a sign change under a CCW exchange, while the 
vortex o n the right does not [c.f. Eqs. (BlOc) and (BlOd)]. Thus, the signs of the numerically-obtained overlaps 


(BlOa)-(BlOd) are consistent with the analytic predictions obtained from the action of the braiding generators (8). 


However, note that the magnitudes of the overlaps are not 0.5, as would be expected for perfectly adiabatic braids. 

We identify two reasons for these deviations from ideal behavior. The first is non-adiabatic loss to bulk modes 
of the lattice, which decreases the magnitude of all overlaps between the initial and final states and results from a 
failure to satisfy adiabaticity (see Appendix E.2 for a detailed discussion of the adiabatic condition). The second 
source of non-ideal behavior is adiabatic loss between vortices, which accounts for the differing magnitudes between 
left and right overlaps of the same braid. The two vortex zero modes form a degeneracy-2 subspace of the total 
Hilbert space of the system, and in principle nonzero Berry’s matrix elements between these modes (or equivalently, 
between the symmetric and antisymmetric superpositions of these modes) are allowed. Such matrix elements 
would lead to a final state with left and right decompositions of different magnitude, consistent with the following 
numerically-computed norms: 


{^cw,l\'iPcw,l) = 0.56, 

(Blla) 

(V’cw.flh/’C'w.i?) = 0.44, 

(Bllb) 

{i’ccw,L\i’ccw,L) = 0.42, 

(Bile) 

{i’CCW,R\i’CCW,R) =0.58. 

(Blld) 


[Note that Eqs. (Blla) and ( |Bllb ) add up to 1, as do Eqs. (Bile) and ( Blld| ), as expected.] We expect the Berry’s 
matrix elements leading to the transfer of intensity between vortices to vanish when the vortices are sufficiently far 
apart. However, this requires the vortices to traverse longer distances, which hinders the adherence to the adiabatic 
condition. In principle this could be compensated by increasing the time T of the braiding process or the system 
size; in practice, this leads to a substantial increase in the runtime of the simulation. Note that neither of these 
sources of non-ideal behavior affects the phase of the computed overlaps, which remain consistent with the overlaps 
predicted by the braiding generators. 


4. Derivation of the braiding generators in Eq. (8) 


In this section, we show that the operators 

_ g7r(SJ . 6o,i+i)/2 ^ 

defined in Eq. (8) of the main text, implement the braiding operations 

Ti bo,i fj = bo,i+i 
Ti bo,i+l T- = —6o,i 

^o,j rj = boj , j + l 


(B12) 

(B13a) 

(B13b) 

(B13c) 


explored above. The third line (B13c) follows because if j ^ i or i 1, th en [5 q f^] = 0, since Ti in (B12) only 

contains bosons with labels i and i 1. Note that one can use Eqs. (B12) and ( |B13 ) to show that the braiding 
generators and tj do not commute for all i and j, thus demonstrating directly that the form a non-Abelian 
representation of the br aid g roup. 

To show lines (B13a) and (B13b), it is convenient to make the change of variables 


A,± = (feo,i ±ifeo,i+i) , 


(B14) 
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writing 




\/2i 


(B15) 


Notice that, as defined, [/5i,±,/dj j_] = 1, and $]^^] = 0. In terms of the new bosonic operators /5i^±, we can 
write 

fj = (/3j,+A,+-/3l,_ ft,-) _ (B;L6) 

Now, ^'’+) /3i,+ = e“®5- /3,,+ and /3,,_ = e+®^ /?*,_ e“® ^ , so that 

n /3i,+ r/ = -i/3i,+ (B17a) 

Ti /?*,- r/ = +01- , (B17b) 


?t 


from which (B13a) and (B13b) follow using (B15). 


Appendix C: Waveguide worldlines for vortex braiding 


In this section, we provide explicit formulas for the x-y positions of the waveguides in the presence of v vortices 
in the Kekule order parameter. The positions of the centers of the waveguides in the hexagonal lattice are given by 


va = mi ai + m2 a2. 


rB = rA-\- si, 


(Cl) 


where mi ^2 are integers, the vectors ai ^2 spanning the triangular lattice are given in Eq. (B2), and the vector Si 
connecting sublattices A and B is defined in Eq. (Bl). 


A Kekule distortion, constant or spatially varying, displaces the equilibrium positions of the waveguides to new 
positions va^b + with 


UA{r) = -jA{r)e 


—irK 


-l^l+CX. 


(C2) 


UB{r) = -7 A(r)e 


—irKj^ 


+ C.C., 


(C3) 


where 7 is a dimensionful constant with unit s of [length]^. The pattern of waveguide displacements in the presence 
of V vortices is obtained by substituting Eq. ( |B 6 [ ) into the above. 

In this work, we have modified slightly the above pattern of distortions in the presence of v vortices, in order to 
mitigate the effect of the position-dependence of the waveguide dispersion on the braiding of the associated 
zero modes. The modified pattern of distortions used to ge nerate Eig. 1(B), which leads to a hopping modulation 
of the form (B5) with an order parameter of the form (B 6 ), is 


UA{r) = - 7 A(r)e 


+ c.c. 


UB{r) = 0 . 


(C4) 


We show in Appendix D that, for small displacements ua,b^ this choice of lattice distortion leads to an onsite 
potential that is uniform for all sites in sublattice A, so that the localized states in the vortex core, which have 
support only in sublattice A, experience an onsite potential that is uniform everywhere. (Eor a vortex with the 
opposite vorticity, whose associated zero mode has support on sublattice B rather than A, one would need to freeze 
sublattice A instead of sublattice B.) 

To write into some bulk medium a photonic lattice that implements an arbitrary braid 5, one simply promotes 
the displ acements by introducing a ^-dependence of the vortex positions Ri = Ri{z) in 

Eq. (|B 6 |). The Ri{z) can be defined piecewise on intervals of length Lb (see Appendix E) for each segment of 


the braid corresponding to the application of a single braiding generator. On each of these intervals, Ri{z) can be 
chosen to be an arbitrary parametric function of z that interpolates between the initial and final positions of the Tth 
vortex during that segment of the braid (albeit sufficiently slowly that the adiabatic condition holds). Examples of 
such piecewise parametric functions are shown in Eig. 2 of the main text, where the strands being braided represent 
the vortex-core positions Ri{z) for i = 1, 2, 3. 
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| 0 ) + | 1 ) ^ | 0 ) + | 1 ) 

a) x/2 ® 


| 0 ) - li’) 

V2 


|o) + li’) 

V2 



|0) + IV^) ^ |O) + I'0) 

b) \/2 \/2 


FIG. 4: Simulating the braiding of photonic zero modes with Majorana zero modes. To reproduce the result of a clockwise 
exchange of two photonic zero modes [panel (a)], it is necessary to perform a clockwise exchange of two pairs of Majorana 
zero modes, followed by a full 27r-braid of one of the pairs [panel (b)]. The two processes coincide upon identifying the states 
I 0) and I 1) in panel (a) with the states | 0) and | ) in panel (b). 


Appendix D: Comparison with Majorana zero modes 


The action ( |B13 ) of the braiding generators (B12) on the zero-mode annihilation operators 6o,i is identical to 
the action of the corresponding braiding generators for Majorana zero modes (c.f., e.g., Ref. [26]). The question 
then naturally arises: how similar to Majorana zero modes are the photonic zero modes studied in this work? It 
turns out, as we show below, that the braiding of these photonic zero modes can be simulated by braiding pairs of 
Majorana zero modes. The converse of this statement is not true — photonic zero modes cannot simulate Majorana 
zero modes. Nevertheless, their braiding is still non-Abelian, as evidenced by the noncommutativity of the braiding 
generators (B12). 

To understand how to simulate the braiding of photonic zero modes with Majoranas, first recall that there is a 
well-defined notion of parity for the photonic zero modes: a single zero mode can be occupied by either an even 
or an odd number of photons, and this even-odd distinction affects the outcomes of braiding processes [c.f. the 
discussion following Eq. (10) in the main text]. For Majorana zero modes, on the other hand, only pairs of zero 
modes have a definite fermion parity. This suggests that it should be possible to emulate the braiding of photonic 
zero modes with pairs of Majoranas. 

To see that this is indeed possible, consider performing the braid described by Eq. (11) in the main text-namely, 
the clockwise exchange of two photonic zero modes, each of which is initially in the state (| 0 ) + | 1 ))/a/ 2 , as shown 
in Fig. I^a). (Here, as before, | 0) and | 1) are single-vortex states with zero and one photon occupying the zero 
mode. One could obtain the same result by replacing the state | 0 ) with a state | even) with an even number of 
photons, and the state | 1) with a state | odd) with an odd number of photons.) To perform the same operation 
with Majoranas, one must start with two pairs of Majoranas, each of which is initially in the state (| 0) + | t/; ))/\/ 2 , 
where | 0) is the two-Majorana state with even fermion parity and | ) is the two-Majorana state with odd fermion 

parity. Then, to simulate the exchange of two photonic vortices, one must exchange the two pairs of Majoranas, 
and then perform a full 27r-braid of one of the pairs, as shown in Fig. Qb). As indicated in the figure, the initial 
and final states of the two processes are identical upon identifying the states | 0 ) and | 1 ) in the photonic system 
with the states | 0) and \^p) in the Majorana system. In other words, we have shown that the set of operations 
that can be performed by braiding photonic zero modes is a subset of the set of operations that can be performed 
by braiding Majorana zero modes. This subset is defined by imposing a constraint whereby Majoranas can only be 
braided in pairs, with each braid being of the type shown in Fig. 

This finding has implications for the utility of photonic zero modes in applications for topological quantum 
computing (TQC). In particular, since the braiding of Majorana zero modes is not universal for TQC, [37| it 
immediately follows that photonic zero modes are also not universal for TQC. Furthermore, photonic zero modes 
can perform only a subset of the set of topologically-protected gates that Majoranas can. Thus, photonic zero 
modes do not offer an advantage relative to Majorana zero modes for the purposes of TQC. However, it is a central 
message of this paper that the control and tunability afforded by photonic platforms presents a marked advantage 
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over solid-state devices in terms of the ability to create, manipulate, and observe topological defects and their 
non-Abelian braiding. 


Appendix E: Justifications for neglecting nonuniversal effects 

In this section we discuss the requirements imposed by dynamical phases, adiabaticity, wavefunction mixing, and 
finite size effects, and use these to estimate the number of braids we can perform in our system. 


1. Dynamical phases 

In our tight-binding Hamiltonian, the sublattice symmetry (SLS) ^pr —"^r+sj ? with r in sublattice 

A, ensures that a zero mode will be at precisely zero energy. However, the presence of a position-dependent 
waveguide dispersion (referred to below as an onsite potential) or next-nearest-neighbor hopping breaks SLS 
and can cause the zero mode energy to depend on both a and the position of the vortex center. This could cause the 
dynamical phase accumulation of the zero mode to depend on the path taken by the vortex, foiling measurement 
of the geometric phase due to braiding. We show that, and estimate the total dynamical phase accumulated in our 
system. 


a. Onsite potential 

We show here that the onsite potential can be made uniform in sublattice A by “freezing” (ie., not displacing) 
the sites on sublattice B. Suppose that k^^^r = fhe absence of displacements of the waveguides from the 

hexagonal lattice sites va and r^. When the waveguides are displaced, one can show that, for r in sublattice A, 


kujA = kuj exp ^ 

i=i 


-UA{r) +UB{r + Sj)) 


- 1 


(El) 


For r in sublattice H, one simply takes A ^ B and Sj 
leading order for small displacements, i.e. \ua,b\ ^ n: 


to —Sj in the above. 


Let us now expand this expression to 


CL 

J = 1 


(E2) 


Noting that — 0 [^•/* E)q. (BI)], we see that k^ 

= 0, as in Eq. (C4), renders k^. 


Ub 


= kuj for 


depends solely on ub in this limit. Therefore, setting 
r in sublattice A, to leading order in the displacements of the 
waveguides. (Likewise, setting ua = 0 renders k^^r = for ^ in sublattice B.) Keeping this in mind, we proceed 
to second order in the expansion of k^^r for UB{r) = 0, finding 


ku 


kr. 


1 


3 

j=l 


1 


- ^{sj-UA{r)y 


(E3) 


If \uA{r)\ is constant as a function of r, then this leading correction is in fact independent of r, i.e. the distortion 
of the waveguide lattice leads to a position-independent renormalized dispersion k^^r ^ k^ Sk^. The leading 
position-dependent correction then appears at third order in |ix^(r)|. A position-independent |n^(r)| can easily be 
achieved if the vortex width Iq is sufficiently small, and if the vortex centers Ri{z) never approach a lattice site to 
within a distance io during braiding. 

It is also worth pointing out that it is possible to choose the magnitudes |n^(r)| of the displacements in such a 
way as t o generate a position-independent k^^r^ so long as we take ub = 0. To see this, note that, if UB{r) = 0, 
Eq. (El) is “ultra-local” in r, as its value at a point r is completely determined by ua{t). It is thus sufficient to 
tune one parameter, namely |n^(r)|, in order to fix k^^r to a fixed value for each r. 
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FIG. 5: Schematic showing the length scales used in Sec. 3b. of this Appendix. Blue circles represent waveguides. 


b. Next-nearest-neighbor hopping 


We now outline how the geometry of the waveguide lattice can be tuned such that next-nearest-neighbor hopping 
t' is stron gly su ppressed relative to nearest-neighbor hopping t. We can estimate t and t' by recalling our Schrodinger 
equation (All) for We expect our Wannier states to decay exponentially with decay length Iq = A/(27r\/Ae), 
where A is the wavelength of light in the material A/(27r) = vjuj. For circular waveguides of radius ro and nearest- 
neighbor spacing a (see Fig. [^, we estimate the hoppings to be 


^-{a- 2 ro)/lo 


(E4) 


and 


l' ^ lg-(V^a-2ro)//o^ 
a 

The waveguide radius ro gives us an additional degree of freedom in choosing t and t'. We see that if 

a - 2ro ~ /o < (V^ - l)a 

holds, one can increase Iq to exponentially suppress the ratio t'/t while only algebraically suppressing t. 


(E5) 


(E6) 


c. Dynamieal phase aeeumulation 

We now estimate the total dynamical phase accumulated during propagation over a vertical distance L. As 
we can exponentially suppress dynamical phases due to next-nearest neighbor hopping, we assume the dominant 
contribution to the dynamical phase is from a varying onsite potential. Consider a disorder potential that at some 
waveguide i takes the form = r]i{z) A, where r]i{z) is dimensionless. The energy shift of the zero mode (in the 
continuum limit) is 


5e{z) = J (fx Sjii{x,z) |'0o(^)P 5 
and the phase shift after propagation past a distance L is 


0 = 



dz 

he 
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The root mean square of the phase shift is 




^xy L iz 

^T — 



where ixy ^ « is the correlation length of the disorder in the plane, iz the correlation length in the paraxial direction, 
and ^ is the horizontal extent of the zero mode. 


2. Adiabatic condition 


Braiding vortices in the vertical direction z must be done slowly enough to satisfy the adiabatic condition. We 
introduce factors of fi and define t = zjc to work in units more familiar from quantum mechanics. The amplitude 
to be in the mode, given that the state starts in the zero mode, is given by 


Cn{t) = - 


f 


dt' 


{n\H\Q) 

Eq — En 




. {n\H\0) 

{Eo - EnY 


^i^{Eo-E„)t 2 


sin 




En)t 


One can estimate the probability Pout(i) = 1 “ |co(i)P to escape the zero mode. Using 1 — |co(t)P = 
we have 


Thus we arrive at 


Pout « {‘JEf ^ 




{0\H\n) {n\H\0) 
{En - Eo)^ 


< 


^^(0|ij|n)(n|ij|0) 


A4 

W 

A4 

w 

A4 






(0|ij2|o) 


{2hf 


Pout < ^ 


In our case, say we move the vortex in space with velocity R, then (0 \H‘^\ 0)c = (0 0)c ^ (A/^)^, 

where ^ is again the horizontal extent of the zero mode. Consider the time it takes to move the vortex a distance 
In the same time, light propagates a distance (in the ^-direction) along the guided mode. Notice that in 
the units that we use above, A is energy, not wavenumber. From our tight-binding Hamiltonian we can estimate 
A ^ hc/^. We can then write 


Pout 


< 

rs-/ 


A4 \rj 

(2hc)2 / e \ 
A4 {lJ 



2 



2 


2 
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or, finally, 


Pout 


< 

rs-/ 



Therefore, the condition for adiabaticity is that the distance along the paraxial direction needed to move the 
vortex a distance ^ in the x^-plane must be large compared to be., 


3. Number of braids 


The total length L needed to do N braids is L ~ 'kN where d is the distance between the vortices. 

(Basically, it takes nd/^ movements of the vortex by a distance and adiabaticity requires that the distance in the 
paraxial direction be for each move.) 

The upper limit on the number of braids is reached when Sc/) ^ tt, hence 


leading to 


TT ^ 6(1) ^ 


T] 

7] 



Vtt N , 


N 


TV T] 


d ilyiz 


Now, the ratio Ud = d/^ ^ 20 guarantees that the zero modes in different vortices do not mix. The ratio 
and we want to keep po ^ 0 . 01 , say. So we can write 


N 


TV Tj 


d (-lyf-Z 


^rj^rid ilyiz ■ 


A few comments are in order. First, notice that, as expected, if there is no disorder (77 —> 0), the number of braids 
is unbounded. Second, even if 77 is non-zero, there is a controlled limit where the number of braids can be made as 
large as desired: one can take the limit of large vortex sizes 

Say one takes a pessimistic estimation using 77 ^ 5%, and ^ ^ 5a. Then, N ^ tt/.S x 5^ ~ 800. For 77 ^ 1 % and 
^ - 10a, N 7 r /.02 x 10^ 160, 000. And for 77 - 1% and ^ - 50a, N 7 r /.02 x 50^ 20, 000, 000. 

One can also estimate the vertical length needed per braid, L/N = Tvud^l-y/po- Taking a ^ 10pm and all other 
estimates as above, we have L/N ^ 3 cm if ^ ^ 5a and L/N ^ 6 cm if ^ ^ 10a. 

Finally, it is noteworthy that the correlation length in the paraxial direction is not necessarily of the order of a, 
it can be much smaller. This is because it comes from fluctuations in writing the waveguides during their growth. 
So iz can be possibly hundredths of times smaller than a, in which case the number of braids estimated above goes 
up by a factor of 100 . 


4. Finite size effects 

In a lattice with edges, each vortex will create both a mode localized at the vortex and a mode localized at 
the edge of the system. If the lattice is small enough, there will be an “energy” splitting between symmetric and 
antisymmetric linear combinations of our zero mode and edge mode. Let us assume the splitting to be proportional 
to t e~~^. To prevent the zero mode from tunneling to the edge (a distance dE away) during the course of braiding, 
we require 




TT 


— e . 

t 


(E7) 
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Appendix F: Beam splitters and interference of coherent states 


In this Appendix, we elaborate on the use of beam splitters to interfere the coherent states of light used in our 
experimental proposal. A beam splitter can be thought of as an optical circuit element that interferes photons 
entering two input channels and outputs the result to two output channels. At the level of single photons, a 50/50 
dielectric beam splitter acts as follows: 


^2 





J_ f h\-ih\ \ _ 
V2 \-i b\+b\) ■ 


(FI) 


The factors of i above arise from 7r/2 phase shifts of the photons in each channel upon reflection within the beam 
splitter. This action on single photons in turn induces an action on coherent states of photons: 


_^ g-(|Aip + |A2p)/2 gAi(SI-iSt)/x/2^A2(-iSl+St)/v/2 |q^ ^ 

Note that the output coherent state is not a superposition of coherent states, but rather a new coherent state with 
a 50/50 admixture of photons from the two channels. 

The interferometer proposed in the main text uses two sets of three beam splitters, one for each vortex core: the 
first splits each incoming laser beam into two equal parts, and the second interferes light from the two branches of 
the interferometer. Let us now follow one of the laser beams through each optical circuit element. We denote the 
initial state of light from one of the three laser beams by |A'-, 0). The effect of the first beam splitter is to split this 
beam in two: 


Ai — i X 2 —i Ai + A 2 


V2 ’ V2 


(F2) 


|A',0) 


^ ZlXi 

V2’ V2 


— I 5 ) 


(F3) 


From here, the light in the coherent state \Xj) enters the waveguide lattice executing the braid BiB 2 ^ while the 
light in the coherent state \Xj) enters the waveguide lattice executing the braid B 2 B 1 . The state of light entering 
the second beamsplitter is then | ± Aj, ±Aj), where the choice of sign on each coherent state depends on the braid 
that was performed on that state. Passing this state through the beamsplitter gives 




zt A j % Xj 

V2 


=pi Aj ± Xj 

V2 


±Aj ^ Xj —i(±Aj ± Aj) 

V2 ’ /• 


(F4) 


Therefore, light exiting one output channel of the beam splitter experiences destructive interference, while light 
exiting the other channel experiences constructive interference. Whether this occurs in the first or second output 
channel depends on the relative sign between Xj and Xj after each beam exits its respective waveguide lattice. 



